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The subsystem compatibility problem, which concerns the question of whether a set of subsystem states are 
compatible with a state of the entire system, has received much study. Here we attack the problem from a new 
angle, utilising the ideas of convexity that have been successfully employed against the separability problem. 
Analogously to an entanglement witness, we introduce the idea of a compatibility witness, and prove a number 
of properties about these objects. We show that the subsystem compatibility problem can be solved numerically 
and efficiently using semidefinite programming, and that the numerical results from this solution can be used 
to extract exact analytic results, an idea which we use to disprove a conjecture about the subsystem problem 
made by Butterley et al. [Found. Phys. 36 83 (2006)]. Finally, we consider how the ideas can be used to tackle 
some important variants of the compatibility problem; in particular, the case of identical particles (known as 
iV-representability in the case of fermions) is considered. 



I. INTRODUCTION 

The fact that quantum mechanical particles can be entan- 
gled with each other leads to a notion of compatibility of quan- 
tum states. Given a finite system of quantum mechanical par- 
ticles, the reduced states of a subset of these particles must 
arise from the partial trace of some density state of the entire 
system, and as a result the subsystem states cannot be arbitrary 
positive operators of trace one. A simple example comes from 
the 'monogamy of entanglement' : If two particles A and B are 
in a maximally entangled state, any other particle C cannot be 
entangled at all with either A or B [lL l2D • 

There are many variants of this problem, but the form which 
we will be most concerned with in this paper is the following 
JH]: Given states of all proper subsystems of a multipartite 
quantum system, what are the necessary and sufficient condi- 
tions for these subsystem states to be compatible with a single 
state of the entire system?. A number of techniques and dif- 
ferent approaches have been used to attack this problem, and 
many partial results have been obtained in recent years. How- 
ever, much of the current works on the problem have concen- 
trated on less general forms of the problem, and also on the 
spectral properties of the involved states. 

Another problem that has been the topic of even greater 
study in quantum mechanics is the separability problem, of 
trying to decide when a multipartite quantum state is entan- 
gled or separable (not entangled). One approach to this prob- 
lem is to use the fact that the set of separable states is a convex 
set, and to use all of the geometric facts that follow from con- 
vexity. This turns out to be a very fruitful approach, and has 
led to the development of many powerful mathematical tools, 
most notably the idea of an entanglement witness, for attack- 
ing the separability problem (see [J] for a review). 

What we aim to do here is to utilise the idea of convexity 
to attack the compatibility problem. We will show that in the 
correct mathematical framework the set of compatible subsys- 
tem states can be viewed as a convex set. We will define the 
notion of a compatibility witness, which can be used to show 
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when a set of subsystem states is not compatible with a state 
for the entire system, and we use these to form the first nec- 
essary and sufficient condition for the compatibility problem. 
However, much as is the case with entanglement witnesses, 
without a full characterisation of compatibility witnesses, in 
practice this theorem is not of any operational use. Nonethe- 
less, we do prove a number of interesting results about the 
nature of compatibility witnesses. 

We will show that the ideas of semidefinite programming 
llH 0] can be used to numerically solve the compatibility 
problem. The inspiration for this lies behind similar applica- 
tions of semidefinite programming to the separability problem 
Mama We will give a brief overview of semidefi- 
nite programming and illustrate how it can be applied to the 
compatibility problem. Furthermore, we will show how ana- 
lytic results can be obtained from the numerical results of our 
solution, and we use it to disprove a conjecture of Butterley et 
al. 031] concerning the necessary and sufficient conditions for 
the compatibility problem involving three qubits (the simplest 
non-trivial case). 

There are other variants of the problem that we will con- 
sider after the main problem. An important situation is 
one where we only have partial knowledge of the subsystem 
states. Another is the case where we have n identical particles 
(bosons or fermions); this problem is of particular importance 
in quantum chemistry in determining the ground state ener- 
gies of a system of identical particles 111 211 . We will be able to 
give solutions to both of these variants using a modification of 
the semidefinite program used to solve the main problem. 



II. THE COMPATIBILITY PROBLEM AND PREVIOUS 
RESULTS 

A. Statement of the problem 

Let us first establish some nomenclature. For any Hilbert 
space K., let B(JC) be the set of bounded operators on JC, and 
let T>(JC) be the space of density states (i.e. positive operators 
of trace 1). 

Let Hi (i = I, ... ,n) be Hilbert spaces with finite di- 
mension di. Define H = Hi <£> Hi <£> . . . <£> H n to be the 
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Hilbert space for an n-partite quantum system, and define 
N = {1, . . . , n}. A set of states on every proper subsystem 
then corresponds to a set of operators pa G B{®ieAHi) for 
every AcJV. 

The compatibility problem asks the following question: 
given pA for all Ad N, what are the necessary and sufficient 
conditions for these states to compatible with a state p G H, 
i.e. 

Pa = T?n\aP (1) 

where Trx represents the partial trace on systems Hi for all 
i G X. We say the set of subsystem states are incompatible if 
no such state can be found. 

Some obvious compatibility conditions arise from the fact 
that two states describing overlapping systems must yield the 
same reduced states when the overlap is traced out. Formally, 
we can state this as follows: if A, B, C C N, A n B = A n 
C = 0, then 

Ttb(paub) = Ttc(pauc)- (2) 

Even in the simplest case of three qubits, there is no known 
solution to this problem, and there is no analytic way of de- 
termining whether an arbitrary set of subsystem states is com- 
patible with a full description of the system (or not). 



B. Summary of current results 

There are very few known necessary criteria for the gen- 
eral form of the compatibility problem as stated above. The 
following theorem is an example of such a condition: 

Theorem 1. Let pa = Tt n \aP (where A C N — 
{1, . . . , n}) be the reduced states of a n-party state p. Then, 
defining 

A = J2 (-1)' A| PA (3) 

ACN 

( where pa here is padded out with identities on the traced out 
systems), then if n is odd, A > 0. If all of the systems are 
qubits, then A < 1 also. 

This was proved in the case of three qubits in Oil, and the 
result in its full generality can be inferred from 111311 . Further- 
more, through a study of the classical analog of the compat- 
ibility problem, it was conjectured in ^ that in the case of 
three qubits, the condition < A < 1 is not only neces- 
sary for subsystem state compatibility, but also sufficient. The 
three party case is also studied in detail in (3, where a num- 
ber of necessary conditions on the spectra of any three party 
state and its two-party reduced states are produced. 

The fact that we only have necessary conditions reflects the 
difficulty in the problem. A set of subsystem states can be 
proved to be incompatible if they violate any of these neces- 
sary conditions. However, the only sufficient conditions we 
have amount to checking whether the subsystem states are 



compatible with a specific state p G B(H), and checking these 
conditions for all possible states is not operationally feasible. 

There do however exist many results in more constrained 
versions of the problem. One question that has been particu- 
larly studied is the following: What are the conditions on the 
spectra of the one-particle reduced states for the full state to 
have a given spectrum? For example, a pure state has a partic- 
ular spectrum, and the necessary and sufficient conditions on 
the spectra of the one-particle states for them to be the reduced 
states of a pure state are derived for any number of qubits by 
Higuchi et al. in IU5I1 . for three qutrits by Higuchi in itlal . and 
for two qubits and a 4-dimensional system by Bravyi in lfl7ll . 
A set of necessary conditions were proved for any number of 
arbitrary systems by Han et al. in IU8I1 . Finally, the most gen- 
eral form of the problem was solved by Klyachko 1 19], using 
methods of symplectic geometry similar to those used by Kly- 
achko to solve a long-standing conjecture by Horn, who con- 
jectured a set of inequalities answering the question 'What are 
the possible spectra of a sum of Hermitian matrices with given 
spectra?'. Related to this is the work of Daftuar and Hayden 
Il20ll . in which the possible spectra of a one-party reduced state 
of a bipartite state are determined. 

As mentioned already, most of this work concentrates on 
establishing relationships between the spectra of the reduced 
and (where one can exist) the total states. In this work we 
intend to move away from this focus and concentrate more on 
the geometric aspects of the problem at hand. 

Other works that establish analytic results about the prob- 
lem include the work of Christandl and Mitichison 12111 . in 
which a connection between this problem and the representa- 
tion theory of the symmetric group is obtained, and the work 
of Jones and Linden [22], in which a system of qubits is con- 
sidered, and it is shown that if we are given all the reduced 
states of more than half the number of qubits, there is almost 
always at most one pure state compatible with these subsys- 
tem states. 

Another important variant of the problem is the case of 
identical particles. This extra condition places more con- 
straints on the density state for the entire system, as pure states 
of bosons (fermions) are symmetric (antisymmetric) under ex- 
change of any two particles. A compatibility problem that has 
received particular attention is one that has become known 
as the A^-representabihty problem - this is the case when we 
have N identical fermions. The necessary conditions on the 
one-party reduced states was solved by Coleman: 

Theorem 2 (Coleman i23"ll . Corollary 4. 3 A). A density matrix 
p is the reduced one-party state of a system of N fermions if 
and only all of its eigenvalues are bounded above by 1/n. 

The necessary and sufficient conditions on the two-party 
reduced state however are still not known, and have been of 
great interest in the field of quantum chemistry, since knowing 
the conditions on the two-party states would allow the calcula- 
tion of the ground state energy of an A-fermion Hamiltonian 
involving 2-body interactions only from the 2-party state and 
an effective 2-body Hamiltonian. This is still a topic of active 
research (see IU2I1 and references therein). 

Finally, there has been some research on the complexity 
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of the compatibility problem. Specifically, it has been shown 
that when we have partial knowledge of the reduced states 
(specifically, a number of reduced states that is polynomial in 
the number of quantum systems, each describing a number 
of systems less than some constant), the corresponding com- 
patibility problem is QMA-complete 1I32I1 and hence also NP- 
hard. It was also shown very recently in l24ll that the specific 
iV-representability problem where we are given the 2-party 
reduced state is also in the complexity class QMA-complete. 
In this paper we will not be concerned with detailed analy- 
sis with respect to complexity classes, but we will discuss the 
computational efficiency and complexity of the algorithms in- 
volved. 



III. CONVEXITY AND THE COMPATIBILITY PROBLEM 

In this section we will use the ideas of convexity on the 
compatibility problem, and introduce the notion of a compat- 
ibility witness, a tool that will allow us to detect sets of sub- 
system states that are not compatible with an overall state. 

A. Compatible states as a convex set 

Let us define Ti^ = ®j<zN.j^kHj (i-e. system k removed 
from the tensor product), and let us define O C <3)f =1 B(Ti^) 
by 

O = { O = (Oi, . . . , O n ) I G B{H [l) ), 0\ = O t }. (4) 
We define an inner product on this space by 

n 

(A,B)=Y / Tr(A l B t ). (5) 

i=l 

This space is a natural space to define the set C C O by 

C = {p = (pi, . . . , p n ) I pj = Trj p for some p E V(H)}. 

(6) 

This is the set of compatible subsystem states, or more pre- 
cisely, the set of vectors of the n states of n — 1 systems ob- 
tained from the partial trace of a full state p (pi corresponds 
to the state on Ti^). The states of a smaller number of sys- 
tems can then be obtained from partially tracing an appropri- 
ate (n — l)-system state. 

The compatibility problem can now be cast in the following 
form: Let cr = (cri, . . . , <r n ) be a vector of the n states of n— 1 
systems. We are then trying to ascertain whether cr lies in the 
set C. This is rather similar to the separability problem - given 
a state p, we are trying to ascertain whether it belongs to the 
set of separable states. This set is convex and compact, and 
we will prove the same fact about C. 

Proposition 1. C is a convex, compact subset of O. 

Proof. We will first prove C is convex. Let p,cr E C. Then 

Pi = Ti\ p, Oi = Tr^ a for some quantum states p, a E B(H). 
Let a E [0,1]. Then 

apt + (1 - a)oi = Tr.((ap + (1 - a) a) (7) 



and since ap + (1 — a)a E T>(H), ap + (1 — a)cr E C. 
Now we will prove C is compact. Define the function R : 

B(H) -f e? =1 S(W») by 

R(p) = {Tr lP ,...,Tr n p). (8) 

By definition, R(D(TL)) = C. Since the trace is a continuous 
function, R is a continuous function. Furthermore, T>(TL) is 
a closed and bounded set, and so by the Heine-Borel theorem 
Il25ll . it is compact. Finally, the continuous image of a compact 
subset is itself compact, and hence R(T>(T-l)) = C is compact. 

□ 



B. Compatibility witnesses 

The idea of a entanglement witness arises from the Hahn- 
Banach theorem. An important geometric fact that follows 
from this theorem is the following B26l 12711 : 

Theorem 3. Let G\ , Ci be disjoint convex closed sets in a 
real Banach space, and let C\ be compact. Then there exists 
a continuous functional f, and a E R such that, for every 
c 1 EC 1 ,c 2 E C 2 , /(ci) < a < /(c 2 ). 

We can use this theorem to prove a result about elements of 
O not in C analogous to the result about entangled states in 
& 

Proposition 2. Let er E O, but cr E" C. Then there ex- 
ists W = (Wi, . . . , W n ) E O such that (W, cr) < 0, and 
(W,p) > 0forall P = (p 1 ,..., Pn ) E C. 

Proof. The one-element set S(cr) = {cr} and C are con- 
vex compact (and hence closed) subspaces of the real Banach 
space O. Since also S(cr) n C = 0, we can apply the above 
theorem. Hence there exists a linear functional g : O — * R 
and a E R such that 

g(cr) <a< g{p) (9) 
for all p E C. By the Riesz-Frechet representation theorem, 

n 

g(X) =J2 r b(ViXi) (10) 

i=l 

where V = (Vi, . . . , V n ) G O. Now, defining W = V 
al/n, where I — (li/d±, . . . , l n /d n ), we obtain that 

n n 

J2MW t a t ) < < ^Tr(V^) (11) 

i=l i=l 

which is the required result. □ 

This result immediately leads to the following theorem: 

Theorem 4. Let cr E O. Then the (n — l)-party states repre- 
sented by cr are compatible with a state of the full system (i.e. 
cr E C) if and only if (W, cr) > for all W E O such that 
(W,p)>0forallpEC. 
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We call any W satisfying (W, p) > for all p £ C a com- 
patibility witness, and we will denote the set of compatibility 
witnesses by W. If (W,p) < 0, the subsystem states p are 
incompatible, and we say W detects the incompatibility in p. 

This is the first example of a necessary and sufficient con- 
dition for the compatibility problem, but (rather like the simi- 
lar theorem for entanglement witnesses) it is not immediately 
operationally useful, because it requires a complete character- 
isation of the entire set of compatibility witnesses. 



C. Properties of compatibility witnesses 

To make the introduction to compatibility witnesses easier, 
we will prove a few results about some properties of compat- 
ibility witnesses. The first shows us how to jump between 
a compatibility witness and a positive operator characterising 
this witness. In all that follows this point, if A £ BiTfy'), and 
B £ B(Hi), then the operator A ® B should be interpreted 
as being a member of B(H) rather than B(7i {i) ) <g> B(Hi) (for 
example, if n = 3, and A = A x <g> A 3 £ B(H {2) ), then A(g)B 
should be interpreted as A\ <g> B ® A 3 £ B{TL)). 

Proposition 3. W is a compatibility witness if and only if the 
operator p(W) £ B(Ti) defined by 



establishing the result. 



□ 



p(W) = ^Wi®l h 

i=l 



(12) 



where li is the identity operator on Hi, is positive. 

Proof. Let p £ C be compatible with a state p. We note that 



J2 MWiPi) = Tr 



J2(Wi®li)p 



Tt(p(W» (13) 



and so the result follows (i.e. if one side is greater than zero, 
so is the other). □ 

However, constructing non-trivial entanglement witnesses 
is not as simple as taking Wi to be positive for all i = 
1, . . . , n, as then Tr(Wj<7j) > for any positive <Xj and so the 
compatibility witness W never detects any incompatibility of 
subsystem states. 

We should note here that the function p is not a one-to-one 
function: there exist O £ O (not necessarily a compatibility 
witness) such that p(0) = 0, and hence W = W + O is 
a compatibility witness such that p(W) = p(W), We will 
return to this point later. 

The structure of p(W) also implies the following result: 

Proposition 4. Let XV £ O ( not necessarily a compatibility 
witness), and let Tj £ B(Hi), with Tr(Tj) = 0. Then, defining 
T = Ti <8> ...®T n , Ti(p(W)T) = 0. 

Proof. Define = ® k ^T k . We see that 

Tr(p(W)T) = It ^(Wi ® li)rj (14) 

n 

= ^ Tr(iy t T«) Tr(Tj) = 0, (15) 
i=i 



This can be reduced to a finite number of conditions by 
taking a basis for the subspace of traceless matrices in each 
Hilbert space Hi, and then forming the tensor product of these 
bases (this yields niLi(^i ~ 1) conditions). 

We can also consider this problem in reverse: given a pos- 
itive operator Z £ B(TL), we may want to ascertain whether 
this operator is of the form Z = p(W). This is all encapsu- 
lated in the following proposition: 

Proposition 5. Let Z e B(TL) be positive, such that 
Tr [Z{Tx <g> . . . <g> T n )] = ifTr{Ti) = 0. Then Z = p(W) 
for some compatibility witness W. 



Proof. Let Bk, m (m — 1, . . . 
operators for B(Hk) satisfying 



, d\) be a basis of Hermitian 



Ifc 1 
Bk, m = -j-; Tr(i?fc. m ) = Sim', Tr(Bk,mBk,n) = -j-5 mn 
dk elk 

(16) 

and define 



B m — <8>fc = i-Bfc,r, 



(17) 



Define D = dim7i = Hi^- The set {B m \ rrii = 
1, . . . , df ; i — 1, . . . , ti } forms a basis for B(H), and so we 
can write 



z = z m B„ 



mei 



= 1, 



,<2?; i = 1, 



where I = {m 
the orthogonality conditions (TToTl). z m 
define the set 



(18) 



i}, and (from 



= DTr(ZB m ). Now 



Z c = {m|m j =l,...,^,i = l,...,n; 
m/j = 1 for some fc = 1, . . . , n} . 



(19) 



For Z satisfying the assumptions in the statement of the 
proposition, we have that z m — Tr(ZB m ) = for any 
m ^ Iq, as these basis elements are a tensor product of trace- 
less operators. Hence Z can be written in the form 



Z — y ' z m B m . 

mefc 



(20) 



When m £ Iq, the basis element B m is of the form Bi ® 1^ 
for some i = 1, . . . , n, and so Z is of the form J^i Wi <8> 1;. 
Hence Z — p(W) for some W £ O, and since Z is positive, 
by Proposition[3] W is a compatibility witness. □ 

Included in the above proof is a method of extracting a suit- 
able W £ W such that Z = p(W). The operators z m B m 
for m e If in ( f20b can be used to construct the operators Wi, 
where i is determined by the component of m that equals 1 . 
Where more than one component equals 1, the operator can 
contribute to any of the appropriate Wi. This freedom repre- 
sents the fact that there is not a unique W £ W satisfying 
Z = p(W), but there are in fact many such witnesses. 

Another method of finding an appropriate W £ W from Z 
is through the following proposition: 
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Proposition 6. Let Z = p(W) for some W £ W. Then 

Z = - E d Zn\a®1a (21) 



where cLa = ILeA d>b and I a = <gn e Ah- 

Proof. Using the principle of inclusion/exclusion, we can ex- 
press ( f20b as follows: 



E(-D fe+1 E 



E 



fc=l ACN,\A\=k \ me I 

\mi — 1 when l£ A 

E ^-T^Tr A (Z)®l A 



\ 



22) 



E 



(_1)|A| 



z 



N\A 



ACN,A=£t> 

which is the required result. 



(23) 
(24) 
□ 



Equation d2Tb allows us to express Z in the form ^i® 
li, and so find a suitable W. For example, in the case n = 3, 
given Z = p(W), a suitable W is given by 



Wi = 
w 2 = 
w 3 = 



1 


li 


E> Z 2 <S 


)1 3 


ll ® ^23 


3d±d 2 d 3 
1 


li 


8> I2 <E 


-Z3 


1 2 ® Zl3 


3d\d 2 d3 
1 


Zi 


d\d 2 

g> 1 2 € 


)1 3 


^12 ® 1 3 


3e?id 2 d3 








ds 



(25) 
(26) 
(27) 



We again see the freedom of choice in W here, as some of the 
terms could appear in more than one of the operators Wi, but 
we have chosen this symmetric expression. This method for 
establishing W from p(W) will be important when we come 
to talk about using semidefinite programming on the compat- 
ibility problem. 



D. Optimal compatibility witnesses 

Given two compatibility witnesses W, W, if W detects 
all the incompatible subsystem states that W does, we can 
discard W and use W instead. More generally, we need a 
notion of an optimal compatibility witness, a witness which 
cannot be replaced by another witness that detects as least as 
many incompatible subsystem states. 

In l28ll . the notion of an optimal entanglement witness is in- 
troduced. As a result, a number of the results and proofs here 
are similar to those given in 12811 . but with some differences 
due to the different structure of this problem. 

Let us take compatibility witnesses W satisfying the nor- 
malisation constraint Tr(p(W)) = 1. We start with some 
definitions: 

• A compatibility witness W is tangential to C if there 
exists p G C such that (W, p) = 0. 



• W is finer than W if (W,p) < => (W',p) < 
(i.e. W detects at least as many incompatible subsys- 
tem states as W), 

• W is an optimal compatibility witness if, apart from 
W itself, no compatibility witness is finer than W. 

The normalisation constraint means that we are fairly com- 
paring witnesses (under the above definition aW is finer than 
W when a^O, even though they clealy detect the same set 
of incompatible subsystem states; in fact, without this normal- 
isation, optimal compatibility witnesses would not exist). 

The first observation of note here is that an optimal witness 
is also tangential: if 



min(W, p) 



^min > 



(28) 



then the witness W defined by W[ = Wi — alj/n is a finer 
compatibility witness, and is also tangential to C at the point 
p € C that minimises the above expression. However, since 



(W,p) =Tr (p(W)p) 



(29) 



(where p is a state compatible with p), then A m ; n is simply the 
minimal eigenvalue of p(W). Hence we have the following: 

Lemma 1. W is tangential to C if and only if p(W) has a 
zero eigenvalue, and W is optimal only if p(W) has a zero 
eigenvalue. 

Given a non-optimal compatibility witness W, we may ask 
how to produce a finer witness W from it. Suppose that the 
minimum eigenvalue of p(W) is Aj^ in > 0, and let P be a 
compatibility witness such that (P, p) > for all p 6 O (not 
just the compatible subsystem states), with maximum eigen- 
value A^ lax . Then, by Proposition|5] the operator 



Z' =p(W)--f^p(P) 



(30) 



defines a compatibility witness W' which is by definition 
finer than W, since (W , p) = (W,p) - ^(P,p) < 
W.p. 

The rest of this section is dedicated to proving the converse 
of this, i.e.: 

Theorem 5. A compatibility witness W is finer than W if 
and only if there exists P such that (P, p) > for all p € O 
and 



W = (1 - e)W' + (P 



(31) 



where e € (0, 1). 



Before we are in a position to prove this theorem, we need 
the following lemma, which we prove here first. 

Lemma 2. Let p € <D,<£ C. Let I = l/D(dili, . . . ,d n l n ), 
the reduced states of the maximally mixed state on Ti. Then 
there exists x € (0, 1) such that p(x) — (1 — x)I + xp is 
compatible with a quantum state p(x) € B(H). 
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Proof. For all compatibility witnesses satisfying 
Tr(p(W)) = 1, we have (W, I) = 1. Define 

c = min(W, p) (32) 
w 

where the minimisation is over all compatibility witnesses W 
satisfying Tr(p(W)) = 1. Since p C, c < 0. Also, 
mm WeW ^ Tr:{p{w))=1 (W , p(x)) = (1 - x) + xc, and be- 
cause of Proposition|4] for x < 1/(1 — c), p(x) is compatible 
with an overall state. □ 



The following lemma is very similar to Lemma 1 in 1 28 
and hence the proof is almost identical, with some small 
changes. We state the proof here however for completeness. 

Lemma 3. For a compatibility witness W, let Dw = {p € 

O | (W, p) < 0}, and for two witnesses W, W, define 



A = min 

crED w 



(W, 



(W,(t) 



(33) 



Suppose W is finer than W, Then the following statements 
are true: 

(i) (W,p) =0=^ {W',p) <0. 

(ii) (W,p) <0^ (W,p) > (W',p). 
(Hi) (W,p) > => X(W,p) > (W',p). 

(iv) A > 1, with equality if and only ifW = W. 

Proof, (i) Suppose (W, p) = 0, and assume (W,p) > 0. 
Take cr S Dw, and define p(x) = xcr + (1 — x)p, so that 
p(x) £ Dw for x £ (0, 1). For small enough x £ (0, 1), we 
have (W, p(x)) > 0, contradicting W being finer than W . 

(ii) Suppose (W,p) < 0. Define a = p + \ {W,p)\I, so 
that (W, cr) = 0. (i) then implies the result. 

(iii) Suppose (W, p) > 0, and take cr £ Dw- De- 
fine p = (W,p)tr + \{W,tr)\p, so that (W,p) = 0. 
Then (i) implies that (W , p) < i.e. \{W, tr)\{W, p) < 
\{W,p)\{W',(t). Dividing both sides by | ( W, <r) | > and 
(W, p) > yields 



(W',p) 
(W,p) 



< 



(W 



(W,a) 



(34) 



Taking the infimum with respect to cr £ Dw on the right hand 
side of this expression then yields the required result. 

(iv) A > 1 follows immediately from (ii). Clearly if W = 
W than A = 1. Conversely, if A = 1, than (i) and (iii) imply 
that (W, p) > (W, p) for all p £ C. Hence if p is a state 
compatible with p, then Tr [(p(W) —p(W'))p] > 0. Let 
Pi = fori = 1, ... , D, with the set {\ipi)}f =1 being 

an orthonormal basis for7i. Since Tr(p(W)) = Tr(p(W')), 
and 

D 

= Tr((p(W)-p(W')) =J2 Tr MW) - p{W')) Pi ] > 0, 

i=l 

(35) 

we must have Tr [(p(W) - p{W'))p t ] = for all i. But the 
orthonormal basis was arbitrary, and this implies that for all 
p £ C compatible with a state p, we have 

(W, p) - (W, p) = Tr [(p{W) - p{W'))p] = 0. (36) 



However, by Lemma|2] for any p £ O, there exists x £ (0, 1) 
such that p(x) = (1 — x)I + xp is compatible with a state 
p(x). Hence (W, p{x)) = (W, p(x)), and hence (W, p) — 
(W',p) i.e. W = W. □ 

Proof of Theorem\5\ First suppose that W = (1 — e)W + 
eP, with (P,p) > for all p £ O. Then (W , p) = 
j^((W,p) -e{P,p)), and hence if (W,p) < 0, then 
(W', p) < also. 

Conversely, take A defined as above. Then Lemma [3 iv) 
tells us that A > 1. If A = 1 then W = W i.e. the theorem 
statement is satisfied with P = 0. If A > 1, then define 



P = 



-(aw - W); 



e = 1 - 



(37) 



1-A V " A 

Then by Lemma Oi)-(iii) above, (P, p) > for all p £ O, 
and by definition, W = (1 — e)W' + eP as required. □ 

An example of a witness P with the above property is one 
where each Pj is positive. But there also exist witnesses P' 
where P[ ^ (since e.g. we can add al to Pi and subtract al 
from Pj , j ^ i to form P' such that (P' , p) = (P, p) > for 
all p £ C. As a result, checking the condition of Theorem[5]is 
difficult, but such witnesses P can easily be used to 'optimise' 
compatibility witnesses as outlined above. 

To conclude this section, we note that there are many more 
results on optimal entanglement witnesses in B28I1 than we 
have presented here. However, the difference in structure in 
compatibility and entanglement witnesses makes it difficult to 
imitate the ideas in [28]. 



IV. SEMIDEFINITE PROGRAMMING AND THE 
COMPATIBILITY PROBLEM 

In this section we outline how the compatibility problem 
can be cast in the form of a semidefinite program, which 
makes it amenable to solution by numerical means. Semidefi- 
nite programming has been used in a number of ways to attack 
the separability problem I^IMI^ fuifTlll. and the ideas utilised 
in those papers (particularly those in 1 10]) act as an inspiration 
for the application here. 

We will start by giving a brief and by no means compre- 
hensive introduction to semidefinite programming. Interested 
readers can refer to ||5j|6j] for more details. We will then out- 
line how semidefinite programming can be used to attack the 
compatibility problem, and we also illustrate how it be used 
to obtain some interesting analytical results, most notably to 
disprove the conjecture in y[]. 



A. A summary of semidefinite programming 

A semidefinite program (SDP) is a minimisation problem, 
which involves the minimisation of a linear function of a vari- 
able x £ W n subject to a linear matrix inequality. Formally, 
the problem is 



Minimise c.x 
subject to P(x) > 



(38) 
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where 

m 

F(x)=F + ^z fc F fe , (39) 
fc=i 

and Fo, Fi, . . . , F m are fixed Hermitian matrices. The min- 
imisation is performed over x £ W n . The set of x for which 
the inequality F(x) > holds is known as the feasible re- 
gion. The SDP is termed feasible if this feasible region is 
non-empty. 

A semidefinite program is an example of a convex opti- 
misation problem, since the feasible region is convex i.e. if 
F(x), F(y) > 0, then for A £ [0, 1], 

F(Ax + (1 - A)y) = AF(x) + (1 - A)F(y) > 0. (40) 

It is this convexity that has allowed fast algorithms to be devel- 
oped to solve semidefinite programs numerically, and details 
about these algorithms can be found in yfl. 

An important structural feature of semidefinite program- 
ming is a feature known as duality. Associated with a problem 
of the form Q8T >. usually known as the primal problem, is a so 
called dual problem: 

Maximise -Tr(F Z) 

subject to Z = Z\ Z > 0, (41) 

Tr(F t Z) = Ci . 

This dual problem is a semidefinite program as it can be cast 
in the form of ( [381 ) [5]. There are a number of relationships 
between these two problems. 

• Let x, Z be feasible solutions to the primal and dual 
problem respectively. Then 

c.x + Tr(FoZ) = Tr(F(x)Z) > (42) 

and so — Tr(F Z) < c.x i.e. any feasible x can be 
used to compute an upper bound on — Ty(FqZ), and 
any feasible Z can be used to compute a lower bound 
on c.x. 

• We say the primal problem is strictly feasible if there 
exits x such that F(x) > (i.e. a strict inequality is 
satisfied). Similarly, the dual problem is strictly feasi- 
ble if there exists Z > satisfying Tr(FiZ) = Cj. If 
either of the two problems are strictly feasible, then the 
optimium values for both problems are equal @]. 

• Finally, if x, Z achieve the optimum values for the pri- 
mal and dual problems respectively, then d42l implies 
that Tr(F(x)Z) = 0. However, since F(x) > and 
Z > 0, F(x)Z = 0. This is known as the complemen- 
tary slackness condition. 

B. Formulating the compatibility problem as an SDP 

We will now discuss how to express the compatibility prob- 
lem as a semidefinite program of the form d38t . The idea is 



this: F(x) will represent the potential compatible state, and 
the inequality F(x) > will correspond to the necessary pos- 
itivity of any state compatible with the (n — l)-party states. 
Within the definition of F(x), the matrix Fo will represent 
the part of any compatible state that is fixed by the reduced 
states. The other matrices represent the freedom in the po- 
tential choice of any compatible state, and as a result do not 
change the partial traces of the potential state. 

Let B m be the basis elements for B(H.) as defined in (fT7l . 
Suppose p £ B(Tt) is a quantum state, with TT k (p) = p^. 
Any p can be written in the form 

P = ^Tp m B m . (43) 

rn 

However, defining Ik = {m | m,k = 1}, and, for m £ Ik, 

B<£ = ®& k B iiTni (44) 
then we can observe that 

Tr k (p) = PmB ( i } (45) 

me/ fc 

i.e. the condition Tvk(p) — p^ fixes p m for m £ I k . Since 
U£ =1 2fc = In (see (fT9il), the n partial trace conditions fix p m 
for m £ Ic 113311 . With this in mind, we define 

B = PmBm (46) 

where the coefficients p m are fixed by the appropriate partial 
trace conditions. 

The remaining basis elements (B m , where m £ I\Ic) are a 
tensor product of traceless matrices (since if m ^ Ic, m k ^ 1 
for all k, and Ti(B k ^ n ) = when m ^ 1), and so yield zero 
under tracing of any system. As a result, the operator 

B(x)=B + ^ x m B m (47) 

mel/Ic 

(where x £ ]R(l- r l - l i °l)) also satisfies the appropriate par- 
tial trace conditions. Therefore, the compatibility problem 
has now been reduced to the problem of finding x such that 
B(x) > 0. 

An equivalent problem is to consider the SDP given by 
Minimise t 

subject to £ + E me j// c Zm£ m + il>0. (48) 

The effect of this problem is to find an operator with the appro- 
priate reduced states with the maximum minimal eigenvalue 
(given by — t). If t > 0, then B(x) can never be positive, and 
the [n — 1) -party states are not compatible with a quantum 
state of the entire system. 

Before we continue, we note that there is a second method 
of calculating the fixed operator Bq: we note that the expres- 
sion for Bq given in d46b is similar to that given in d20b . and so 
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the techniques used in the proof of Proposition|6]can be used 
to deduce that 



Bo 



ACN,A^H) 



where ps (B C N) is the subsystem state for the systems 
labelled by the elements of B obtained by partial traces of 
appropriate [n — l)-party states. 



C. The dual problem and compatibility witnesses 
The dual problem for the compatibility problem is given by 

(50) 



Maximise -Tr(B Z) 
subject to Z = Z\Z > 0, 

Tr(Z) = 1, 

Tr{B m Z) = for all m € I \ I c - 



However, Z = p(W) for some compatibility witness W: 
Since I\Ic = { m I m k = 2, . . . , d 2 ,, fc — 1, . . . , n}, and 
Bk,2, ■ ■ ■ , Bfc^p forms a basis for the space of traceless oper- 
ators on Hk, Proposition 5 implies that Z = p(W) for some 

W eW. ' 

The upshot of this is that if the optimum value of the dual 
problem is positive, the Z which realises this maximum is 
such that Z = p(W), and 



> Tt(BqZ) = Tr 



B ( W i ® !< 



vfe=l 



fc=i 
= (W,p) 



(51) 

(52) 
(53) 



i.e. W is a compatibility witness that detects the fact that 
the (n — l)-party reduced states pi are not compatible with a 
quantum state of the overall state ll34ll . 



D. Complexity 

The SDP (l38l involves optimising over the space R m , 
and dealing with the positivity of N x N Hermitian ma- 
trices. Numerical SDP solvers are known to have com- 
plexity 0{m 2 N b / 2 ) ((^(iV 1 / 2 ) iterations, each taking time 
0(m 2 N 2 ) @]) 01. 

In our case, if we define D — JT, <h, then 

• The number of variables is equal to |/| — \Ic\ + 1 = 

nr=i(^-i) + i-o(^ 2 ) ; 

• The dimension of the matrices equals D. 

and so the complexity of the compatibility SDP is O(Z) 13 / 2 ). 
This is polynomial in the dimension of every individual sys- 
tem, but is exponential in the number of parties i.e. if dj = d, 



then D — d n , leading to exponential complexity in n. How- 
ever, the complexity of the algorithm is still polynomial in 
the problem size, since specifying the reduced state which ex- 
cludes system k involves specifying a D/dk x D/dh Hermi- 
tian matrix i.e. (D/dk) 2 parameters. So in total we have a 
problem size of 



O D- 



fe=i 



0(D 2 ) 



(54) 



hence the time complexity is bounded above by a polynomial 
in s. Hence the problem lies in the complexity class P. The 
exponential complexity in n merely reflects the fact that we 
must supply an exponentially increasing amount of initial data 
in the form of the reduced states. 

Note that this does not contradict the results of ll29ll dis- 
cussed earlier, as the compatibility problem considered there 
involves having only a number of reduced states that is poly- 
nomial in the number of quantum systems n, and each reduced 
state describes the state of a number of systems below some 
fixed constant. Hence the problem size here is only a polyno- 
mial in n, leading to very different conclusions regarding the 
complexity of the problem. We will return to this issue later. 



E. New insights from the compatibility SDP 

We have produced computer code for the compatibility 
SDP using MATLAB and a freely available SDP numerical 
solver known as SEDUMI lf30h . In much the same way as 
was achieved when semidefinite programming was applied to 
the separability problem, analytical results can be gained by 
inspection of numerical results due to the accuracy numerical 
SDP solvers work to. 

A simple example of this comes from considering the fol- 
lowing two-qubit subsystem states of a three party system: 



P 



(i) 



,(2) 



-,(3) 



1 



|*+)<*+|; |*+) = — (|00> + |11» 

(55) 

which are not compatible with an overall state (if qubits 1,2 
are in a pure entangled state, then qubit 3 must be in a pure 
state, which is not the case here). The primal SDP confirms 
this incompatibility, and the dual SDP can be used to deter- 
mine a numerical compatibility witness. From this result, an 
analytic compatibility witness can be derived that detects this 
incompatibility. Let 



^ = ^(1^1)^11 + 1^X^1) 



(56) 



where 



IV»i) = 2(|000)-|011)-|101)-|110)), (57) 

|V>2> = ~(|001) + |010) + |100)-|111>)- (58) 

This can easily be checked to be a compatibility witness by 
using Proposition [6] Due to the symmetry of Z, all the three 
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two-party reduced states of Z are equal: 



1 



z u = z 13 = z 23 = - (|*_X*_| + |*+)<*+|) ; (59) 



Let a — j (|1)(1| ® 1 ® 1). Let us define the following 
two-qubit subsystem states of a three qubit system: 

pW(p) = (l-p)Tr fc <7+p|* + )(f + | (63) 



V2 



(|00)-|11)),|* H 



V2 



(|01> + |10)). 



(60) 

and all the one-party reduced states of Z are equal to 1/2. 
Hence a compatibility witness W corresponding to Z is given 
by 



W 1 = W 2 = W 3 = + I (|*_)(* 



■ l$ + )<*+l) 

(61) 

1/4, detecting the in- 



which leads to ^Li Tr(p^Wi) = 
compatibility as required. 

The compatibility SDP can be used to help us produce a set 
of 2-qubit reduced states that disproves the conjecture in fl. 
For the case of 3 qubits, the conjecture states that ,012, pi 3 , P23 
are the 2-qubit reduced states of a three-qubit state p if and 
only if 

< A(p) = 1 - pi - p 2 - p 3 + P12 + p%3 + P23 < 1 (62) 

where the reduced states are padded out with identities as 
stated in TheoremQ] and p is the 3-vector of 2-party reduced 
states. 



where p 6 [0, 1]. This set of subsystem states is clearly com- 
patible with an overall state when p = 0, and incompatible 
with any overall state whenp = 1. We can calculate A(p) for 
these reduced states to be 

A(a) = (1 - p)l/4 + 3p/2 (|^iX0i| + ]^2><^a|) 

- p/2 (| + \ih){ih\) (64) 



where 



n = 



\<k 



(3|000) + |011) + |101) + |110)), (65) 
(|001) + |010) + |100) + 3|111)), (66) 



and hence < (ip\A\ip) < 1 whenp < 1/3. If the conjecture 
is true, then when p < 1/3, there should be an overall state 
compatible with the above subsystem states. However, the 
primal SDP can be used to show this is not true, and produce 
a numerical compatibility witness illustrating this fact. From 
this, an analytical compatibility witness can be obtained. For 
the case p = 1/4, it comes from the operator 



Z 





/ 2959 








-102 





-1715 


-1715 


° \ 







24865 


1005 





766 








-1715 







1005 


24865 





766 








-1715 


1 


-102 








45033 





766 


766 





100000 





766 


766 





46 








-102 




-1715 








766 





1006 


1005 







-1715 








766 





1005 


1006 







V 


-1715 


-1715 





-102 








228 / 



(67) 



which can be verified to be positive, to be of the form Z = 
p(W) through Proposition 6, and to witness the incompatibil- 
ity of the above subsystem states for p = 1/4. This example 
hence disproves the conjecture given in JH. This witness can 
also be seen to be an example of the optimisation process we 
saw in section HIID I Defining 

Z' = 10" 6 diag(0, 0, 0, 0, 2913, 23859, 23859, 44805) + al 

(68) 

and considering the decomposition Z = (Z + Z') — Z' , then 
it can be shown that 

(a) Z + Z' = A(Z + Z'); 

(b) Z'=p(P); 

(c) For large enough a, (P, p) > for all p € O. 



i.e. Z is an optimisation of A(Z + Z'), and can detect incom- 
patibility that the condition A(p) > cannot. 



SOLUTIONS OF VARIANTS OF THE COMPATIBILITY 
PROBLEM 



This new technique of using semidefinite programming to 
numerically solve the compatibility problem can be applied 
to important variants of the compatibility problem we have 
considered. We discuss these applications in this section. 
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A. Partial knowledge of reduced states 

A simple modification of the problem is to consider a situ- 
ation where we are only given a partial subset of the reduced 
states. The problem (f48b can easily be modified to solve this 
problem: The number of basis operators B m that appear in 
the fixed operator Bq is reduced according to which reduced 
states are known, and the coefficients of the remaining basis 
operators are then determined by the SDP solver. The dual 
problem to this program will still produce a compatibility wit- 
ness Z = p(W), but with some extra conditions of the form 
Tr(ZB m ) = for the the basis operators that are no longer 
fixed in the operator Bq. 

The complexity of the problem depends on the problem size 
i.e. the number of known reduced states, and the number of 
systems described by these states. Let us consider the case 
discussed in j29tl (to reiterate: we have a number of reduced 
states p(n), a polynomial in n, each describing a number of 
systems less than some constant fc). In this case, the prob- 
lem size is 0(p(n)), whereas the dimension of the matrices 
involved is D, and there are 0(D 2 ) free variables (since fix- 
ing polynomially many reduced states describing at most fc 
systems can only fix polynomially many coefficients in Bq). 
Hence the complexity is again 0(D 13 / 2 ), but this is no longer 
polynomial in the problem size. 

This is simply one complexity scenario we could imagine; 
we could place other conditions on the number of reduced 
states we know, and how many states they can describe. As 
complexity is not our main focus here, we leave these ques- 
tions for further investigation. 



B. The bosonic and fermionic compatibility problems 

If a density state p represents the state of either n bosons or 
n fermions, there are some additional constraints that much 
be placed on p. A pure state \ip) e TL representing n bosons 
must be totally symmetric under particle exchanges (or swaps 
between Hilbert spaces); whereas any pure state of n fermions 
should be totally antisymmetric under particle exchanges. 

Since we are dealing with identical particles, the n single 
particle Hilbert spaces are all identical; let their dimension 
be d. Also, let -Bfc,m be equal for all k to B m . The total 
(anti)symmetry means that all of the fc-party reduced states 
will be identical; as a result we will refer to any fc-party re- 
duced state as the fc-particle reduced state. 

The space of totally symmetric pure states has an orthonor- 
mal basis {\S x )} 1 < Xl <...< Xrl < d , where 

/ ni(x)!...n d (x)!\ 1/2 ^ 
\S*)={ - } I 2^ \x p{1) )...\x p(n) ) 

(69) 

(where (x) is defined to be the number of instances of fc 
within x). The projection on this subspace is given by 

Ps= E \S X )(S X \ (70) 

xi < . ..<x n 



A density state p is then totally symmetric if and only if 
PspPs — P- However, since p is positive and of trace one, 
this relation is satisfied if and only if Tr(PspPs) — Tr(p) = 
1. Finally, since P| = P$, and by the cyclic property of 
the trace function, we have that p is totally symmetric if and 
only if Ti(pPs) = 1. Noting that Tr(P s pP S ) < Tr(p) = 1 
(because Ps is a projection), we can change the condition to 
Tr (pPs) > 1- Taking all this into account, we can hence 
modify problem (l48b to accommodate the totally symmetric 
condition: 

Minimise t 

subject to Bq + J2mei\i c Im5m + *1 > 

-1 + Tr(B P s ) + Eme/Vc Tr ( S mPs) > 

(71) 

(note that two linear matrix inequalities B 1 (x) > 0, B 2 (yi) > 
can be combined into one of the form B\ (x) © B 2 (x) > 0, 
so this really is a semidefinite program). 

We can produce a similar SDP for the fermionic case. The 
space of totally antisymmetric pure states is spanned by an 
orthonormal basis {\A x )}i< Xl< ... <Xri <:d, where 

= E £ (p)Kw) ■ ■ ■ K(n)) ( 7 2) 

' pes„ 

where £(p) denotes the sign (parity) of the permutation p € 
S n . The projection onto the antisymmetric subspace is hence 
given by 

p A = J2 \ A *)( A *\ < 73 > 

Xx<...<X n 

and similarly the semidefinite program for the fermionic com- 
patibility problem is given by 

Minimise t 

subject to B + J2mei\i c I "> S "> + tl > 

-1 + Ty(B P a ) + E m ei\io x ™^( B ™ p A) > 

(74) 

We can reduce number of free variables in both of these 
problems by noticing that for both the bosonic and fermionic 
case, the density matrix of n particles is symmetric under 
any particle exchange (in the antisymmetric case, |A X ) is to- 
tally antisymmetric, but any operator |^4 x )(A y | is symmet- 
ric because any factors of —1 cancel out). Hence the den- 
sity matrices only depend on symmetrised basis operators 
{- B m}i<m 1 <...< m „<d, where 

B m= E \<S> B rn p(k) ) ■ (75) 

pgs„ \k=i J 

Further to this, the two semidefinite programs can be modified 
again to deal with cases where we have only partial knowl- 
edge of the subsystem states, by adjusting the inequalities in 
the same way as outlined in section [VA] This would allow 
us for example to deal with the iV-representability problem 
discussed in section HTbI 
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C. Determination of ground state energies of identical particle 

systems 

We can further modify the semidefinite programs fi7T[ and 
( l74i i to create a program for determining the ground state 
energy of a system of n identical particles with a Hamil- 
tonian that consists only of two-body interactions. The 
(anti)symmetry of any state of these particles implies that we 
can replace the 7i-body Hamiltonian with an effective two- 
body Hamiltonian H^ 2 \ The 2-particle reduced state p^ 
of a state p can then be used to evaluate the energy E — 
Tr(i?( 2 ) pW) (see e.g. JH for more details). 

We can use semidefinite programming to minimise E, sub- 
ject to constraints that ensures is a valid 2-particle den- 
sity state. The following program calculates E for a set of n 
bosons: 

Minimise Tr(H^p^) 
subject to p (2) = Tr 3 N p (76) 
p>0;Tr(pP s )>l 

(by writing p = Y,i< mi <...< mn <d x ™ B L *e above prob- 
lem can be cast in a semidefinite program of the form of d38l)). 
For a set of n fermions, we simply replace P$ by Pa- How- 
ever, the complexity of this problem will again be exponen- 
tial in n, and so for large numbers of particle this solution 
will be inefficient. Indeed, recent results showing that the 
TV-representability problem is QMA-complete 12411 illustrates 
that an efficient solution to this problem (one with running 
time polynomial in n) is highly unlikely to exist. 

D. Other variants 

The key fact that allows us to apply the ideas of semidefinite 
programming to these variants of the subsystem problem is the 
fact that both the set of (anti)symmetric states and the set of 
any fixed subset of the subsystem states compatible with an 
overall state are both convex sets. A variant of this problem 
that would be difficult to solve using this algorithm would be 
the problem where we wish to find a compatible state with a 
particular fixed spectrum, because the set of operators with a 
fixed spectrum is not convex. A different method would have 
to be used in this case. 



VI. FURTHER IDEAS 

In this paper we have illustrated how we can numerically 
solve the compatibility problem using semidefinite program- 
ming, and the fast algorithms available to solve semidefinite 
programs make such a task operationally feasible l36ll . With 
only slight modifications, the SDP can deal with a situation 
where we only have partial knowledge of the proper subsys- 
tem states, or where we are dealing with a set of identical 
bosonic or fermionic particles. 

The introduction of compatibility witnesses in this paper 
marks a move away from the usual methods used to attack the 



compatibility problem, and as a result there is possibly more 
that can be said about the problem. We will discuss some of 
those possibilities here. 

Having disproved the conjecture from [3], we are left then 
with what further necessary conditions we can find for the 
compatibility problem, and what (finite or infinite) subset of 
these form a sufficient condition for the problem. One possi- 
ble source of these could come from a study of positive maps: 
We can view the operator A as a positive map on the space 
B(TL) that only depends on the reduced states pa- So, find- 
ing other positive maps only depending on the reduced states 
could give us more necessary conditions. Finding whether 
these conditions could be sufficient however would be diffi- 
cult. 

Another possible source of ideas is a study of the semidefi- 
nite program itself. The problem of minimising the maximum 
eigenvalue of a symmetric matrix that is an argument of an 
affine parameter has been well studied (see references in (Hi), 
and these ideas could lead to new insights about this problem. 

One oddity about this problem is the fact that many of the 
useful results obtained so far only apply when we have an odd 
number of systems. We can obtain conditions for the even 
case by treating two individual systems as one system, but 
this obviously does not reflect the full generality of the prob- 
lem. Despite some effort, we have not been able to make any 
further progress in this case. 

Finally, we presented a number of properties of compatibil- 
ity witnesses, but we have not given a complete characterisa- 
tion of these witnesses, and knowing more about their struc- 
ture would be very helpful indeed. However, we will make a 
further observation here. Let us continue to explore the case 
where we have an odd number of systems (n is odd). For 
Z G B(H), let 

A(Z)=J2(- 1 ) IAIz a®1n\a (77) 

AcN 

(This is equal to the expression defined in Theorem Q] when 
an overall state exists). Then if Z > 0, A(Z) = p(W): Let 
T = Ti <g> . . . ® T n , with Ti G B(Hi), Tr(T,») = 0. Then 

Tr(TA(Z)) = M)'^ MZaTa) Tr(T N \ A ) = 

ACN 

(78) 

and since A is a positive map [13], Proposition [5]imples that 
A(Z) = p(W). However, this is does not form the entire set 
of compatibility witnesses: For p ^ C, but A(p) > 0, then, 
for all Z > 0, 

Tr(B A(Z)) = Tr(A(B )Z) (79) 
= Tr(A(p)Z) > (80) 

where the first line follows from the fact the adjoint map of A 
is itself, and the second line follows since the partial traces of 
Bq coincide with the subsystem states in p. All of this implies 
the following theorem: 

Theorem 6. Ifp^C, but A{p) > 0, then for all compatibil- 
ity witnesses W satisfying p(W) = A(Z) for some Z > 0, 
(W,p)>0. 
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Again using our analogy with the separability problem, we 
can think of sets of subsystem states that are incompatible 
with a full system state, but satisfy A(p) > as analogous to 
entangled states with a positive partial transpose, in the sense 
that it is difficult to explicitly find entanglement witnesses that 
detect these states. The compatibility witness given at the end 
of the last section is an example of a witness that does detect 
incompatible subsystems states satisfying A(p) > 0, and fur- 
thermore it comes from an enhancement of a witness of the 
form A(Z). 

There is great potential for powerful compatibility wit- 
nesses to be constructed using this witness optimisation tech- 
nique. Exploring this avenue further and looking for other 



methods to construct compatibility witnesses not of the form 
A(Z) could well be a good step into producing new necessary 
criteria for the problem. 
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